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Abstract 

Epigenetic regulation of gene expression has been shown to change over time and may be associated with environmental 
exposures in common complex traits. Age-related hearing impairment is a complex disorder, known to be heritable, with 
heritability estimates of 57-70%. Epigenetic regulation might explain the observed difference in age of onset and 
magnitude of hearing impairment with age. Epigenetic epidemiology studies using unrelated samples can be limited in 
their ability to detect small effects, and recent epigenetic findings in twins underscore the power of this well matched study 
design. We investigated the association between venous blood DNA methylation epigenome-wide and hearing ability. 
Pure-tone audiometry (PTA) and lllumina HumanMethylation array data were obtained from female twin volunteers enrolled 
in the TwinsUK register. Two study groups were explored: first, an epigenome-wide association scan (EWAS) was performed 
in a discovery sample (n = 115 subjects, age range: 47-83 years, lllumina 27 k array), then replication of the top ten 
associated probes from the discovery EWAS was attempted in a second unrelated sample (n = 203, age range: 41-86 years, 
lllumina 450 k array). Finally, a set of monozygotic (MZ) twin pairs (n = 21 pairs) within the discovery sample (lllumina 27 k 
array) was investigated in more detail in an MZ discordance analysis. Hearing ability was strongly associated with DNA 
methylation levels in the promoter regions of several genes, including TCF25 (cgOl 161216, p = 6.6x10~ 6 ), FGFR1 
(cg15791248, p = 5.7x10" 5 ) and POLE (cg18877514, p = 6.3x10" 5 ). Replication of these results in a second sample 
confirmed the presence of differential methylation at TCF25 (p(replication) = 6x10~ 5 ) and POLE (p(replication) = 0.016). In 
the MZ discordance analysis, twins' intrapair difference in hearing ability correlated with DNA methylation differences at 
ACP6 (cgOl 377755, r= -0.75, p=1.2x10~ 4 ) and MEF2D (cg081 56349, r=-0.75, p=1.4x10~ 4 ). Examination of gene 
expression in skin, suggests an influence of differential methylation on expression, which may account for the variation in 
hearing ability with age. 
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Introduction 

The term epigenetics [1] refers to the regulation of gene 
expression primarily by DNA methylation and changes to DNA 
folding. Epigenetics plays an important role in gene expression 
regulation and cell differentiation in the developing organism 
[2,3] . While the genetic code is fixed, epigenetic changes may be 
dynamic and have been shown to change during a lifetime [4]. 
DNA methylation is one of the most commonly studied epigenetic 
changes and involves the addition of a methyl-group to the 5th 
carbon molecule of a cytosine base, generating 5-methyl-cytosine. 
This stable modification occurs primarily at the CpG dinucleotide, 
but has also been detected at CpH sides, where H can stand for C, 
A or T. Each diploid human genome contains on average 10 



cytosines, of which about 1 0 are combined with guanine as CpG 
dinucleotide [5]. CpG dinucleotides often cluster in CpG-islands 
in the promoter region of genes. The majority of CpG islands in 
promoters are unmethylated. DNA methyl-transferases are 
responsible for de novo methylation of DNA [6] and facilitating 
stable transmission of epigenetic marks during cell division [7] . 

Epigenetic changes can be influenced by both environmental 
exposure [4,8] and genetic variation [9]. Therefore, samples for 
studying the association of methylation with any given trait should 
ideally be matched for genetic and environmental variation. This 
could best be achieved by using family data or monozygotic twins, 
which are assumed to be genetically identical and well matched for 
environmental exposures [10]. Family and twin studies have 
identified differentially methylated regions associated with age [11] 
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and multiple complex traits [4,12] and have further been used to 
estimate rates of heritability in DNA methylation [8,1 1]. Current 
advances in technology allow for high-resolution screening of 
DNA methylation profiles across the genome. Multiple platforms 
exists, but to date the majority of studies have successfully used the 
Illumina Infinium HumanMethylation 27 k and HumanMethyla- 
tion 450 k Bead Chips to assay genome-wide DNA methylation 
profiles across individuals [11,13]. 

Age-related hearing impairment (ARHI) is a common complex 
trait affecting 46% of the population over the age of 48 [14]. 
Epigenetic changes in the ageing ear have been proposed to 
account for age-related changes to hearing ability and syndromic 
forms of hearing loss [15,16]. Changes in DNA methylation are 
influenced by environmental exposure and could therefore provide 
the essential link between the environment and changes in gene 
expression. Furthermore, epigenetic changes with age could 
explain how a previously healthy individual develops hearing loss 
with age. Several forms of syndromic hearing loss, such as Rett 
and Stickler syndrome, have been associated with epigenetic 
change [17,18,19]. Here, an epigenome-wide association study 
(EWAS) of hearing ability was performed, the first EWAS of 
ARHI to our knowledge. This research aimed to determine 
significant associations of differentially methylated regions with 
hearing ability in subjects from the TwinsUK cohort. The most 
significandy associated CpG sites were replicated in an indepen- 
dent sample and gene expression profiles were investigated at the 
genes identified. 

Materials and Methods 

Ethics statement 

The study was approved by the National Research Ethics service 
London-Westminster (REC reference number: 07/H0802/84). 
Fully informed written consent was obtained from all participants 
prior to study conduction. All research described was conducted 
according to the rules described in the Declaration of Helsinki. 

Subjects 

Hearing data in form of air-conduction PTA was collected from 
participants of the TwinsUK cohort between 2009 and 2013. 
Hearing thresholds were determined at frequencies 0.125-8 kHz 
for each ear according to the recommendations of the British 
Society of Audiology [20] . Pure-tone audiometry information was 
summarised by principal component analysis (PCA) [23]. All 
participants completed a questionnaire covering exposure to 
environmental risk factors for ARHI and previous ear diseases. 
Subjects reporting a family history of hereditary hearing loss or 
signs of conductive hearing loss were excluded from the analysis. 

DNA methylation profiles 

Whole blood samples for DNA methylation screening were 
profiled using two different DNA methylation assays, the Infinium 
HumanMethylation 27 k BeadChip (26,690 CpG sites) and the 
Infinium HumanMethylation 450 k BeadChip Kit. In both arrays, 
the DNA methylation level at a specific CpG site is expressed as 
the fi value, which represents the ratio of the methylated probe 
signal over the methylated and unmethylated probe signals. The |3 
score ranges from 0 to 1, where 0 indicates absence of methylation 
and 1 represents a fully methylated CpG site. 

The Illumina Infinium HumanMethylation 27 k array measures 
methylation at 27,578 CpG sites, covering 14,495 genes. This 
array covers primarily CpG sites located in promoter regions of 
genes with on average two CpG sites per consensus coding 
sequence and three to twenty assays per cancer gene [21]. The 



Illumina Infinium HumanMethylation450 k array covers 485,577 
methylation sites in 99% of RefSeq genes (21,231 genes) with an 
average of 17.2 CpG sites per gene region. [22]. 

To identify potential confounders of the Illumina Infinium 
HumanMethylation 27 k array, principal component analysis 
(PCA) was performed using the normalised DNA methylation 
values. The first five principal components resulting from this 
analysis were correlated with following covariates: chronological 
age, methylation chip and position of sample on the chip. Both 
methylation chip and position of sample on the chip were 
significandy correlated with the first two principal components 
from this analysis and therefore included as fixed effects in further 
analysis [1 1]. The same procedure was performed for the Illumina 
Infinium HumanMethylation 450 k array with covariates age, 
chip, position of sample on the chip and bisulfite converted DNA 
concentration levels. Chip, position on the chip and bisulfite 
converted DNA concentration levels were significantly associated 
with the first 3 principal components and were therefore included 
as fixed effects in the linear mixed effects models. 

Epigenome-wide association study 

The discovery EWAS of hearing was performed in 1 1 5 adult 
female subjects with available PTA data and Illumina Infinium 
HumanMethylation 27 k profiles [11]. The DNA methylation 
profiles used in the EWAS were obtained from 26,690 DNA 
methylation probes, which mapped uniquely to the human 
genome (hg 1 8) [9] . After further exclusion of probes mapping to 
the X-chromosome and probes with missing data, 24,641 
autosomal probes remained for the EWAS [11]. The Illumina 
Infinium HumanMethylation 27 k profiles have been published 
previously [11]. DNA methylation was transformed to a standard 
normal distribution per probe using quantile normalisation. To 
determine the association between hearing ability and DNA 
methylation a linear mixed effect model was applied. DNA 
methylation levels at each CpG site were regressed against hearing 
ability (PCI), with adjustment for age, methylation chip, order of 
samples on the chip and twin relatedness. To exclude associations 
with DNA methylation due to covariates other than hearing, the 
full model was compared to a null model, excluding hearing as a 
predictor variable. The null and the full model were compared for 
model fit in an analysis of variance (ANOVA). Only associations 
where the full model fitted the data significantly better (p<0.05) 
than the null model were reported. For each significandy 
associated probe, the effect size (beta), standard error of effect 
(se) and the p-value from the analysis of variance comparing full 
and null model were reported. To confirm that the regions of 
association were not age-dependent differentially methylated 
regions (age DMRs), models including and excluding age as a 
fixed effect were compared. Furthermore, associated probes were 
checked against previously reported age DMRs [11]. In addition, 
to exclude an underlying association between genetic (rather then 
epigenetic) variation and PCI, genetic variants in the DMR 
genomic loci were tested for association with PCI, in a PCI 
genome-wide association scan from the TwinsUK cohort 
(n= 1028). To adjust for multiple testing in the EWAS initially a 
Bonferroni corrected significance threshold assuming 24,641 
independent tests (p = 0.05/24641 = 2.03 x 10~ 6 ) was assumed 
epigenome-wide significant. Furthermore, since the Illumina 
Infinium HumanMethylation 27 k array contains on average 2 
probes per promoter and high levels of co-methylation between 
nearby probes have previously been reported [9], we also 
considered 2 additional Bonferroni corrected thresholds: a 
genome-wide significant threshold correcting for 14,495 indepen- 
dent genes (p = 0.05/14495 = 3.45 xlO" 6 ) and a genome-wide 
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suggestive threshold correcting for 14,495 independent genes 
(p = 0. 1/14495 = 6.90 xlO" 6 ). 

Replication study 

The replication sample consisted of 203 females from the 
TwinsUK registry. For the replication study only the 10 probes 
most highly associated in the discovery EWAS were investigated, 
while the remaining 485567 probes from the 450 k array were 
neglected. The 1 0 selected probes were examined for replication in 
the second sample using a linear mixed effect model. DNA 
methylation was transformed to standard normal per probe using 
a quantile normalisation. DNA methylation at each CpG site was 
regressed against hearing ability (PCI) with adjustment for age, 
methylation chip, order of samples on the chip, bisulfite 
conversion levels and twin relatedness. To exclude association 
with DNA methylation due to covariates, the full model was 
compared to a null model, in which hearing was excluded as a 
predictor variable. The null and the full models were compared for 
model fit in an analysis of variance. For each of the 10 probes, the 
effect size (beta), standard error of effect (se) and the p-value from 
the analysis of variance comparing full and null model were 
reported. Replication of association was considered if association 
was in the same direction and nominally significant (p^O.05). To 
confirm that replicating probes were not age-dependent DMRs, 
models including and excluding age as a fixed effect were 
compared. To determine the significance and effect of joint 
association signals in the discovery (27 k) and replication (450 k) 
samples, a meta-analysis was conducted for the ten most highly 
associated probes using METAL [23] based on the inverse- 
variance option. 

DMR validation using methylated DNA 
immunoprecipitation sequencing (MeDIPseq) 

To further validate the findings from the EWAS (27 k) and 
replication study (450 k) using an alternative technique, the top 
ranked DMR was also explored using methylated DNA immuno- 
precipitation followed by high throughput sequencing (MeDIPseq) 
data. The MeDIPseq validation sample consisted of 46 unrelated 
healthy females with PTA scores and previously published 
MeDIPseq profiles [24] . MeDIPseq DNA methylation levels were 
generated and quantified as previously described [24], and relative 
methylation scores in a 1 kb region on chr 16 (chrl6: 88466501- 
88467500 on hg 18) overlapping probe cgOl 161216 (chr 16: 
88466949 on hg 18) were explored for association with PTA. A 
linear fixed effect model was applied, where the DNA methylation 
signal at the locus surrounding the chromosomal position of probe 
cgO 1161216 was regressed on hearing ability (PCI), adjusted for 
age. To exclude an association of DNA methylation with age, the 
full model was compared to a null model, excluding hearing as a 
predictor variable. The null and full models were compared for 
model fit using analysis of variance (ANOVA). 

Whole blood cell subtype heterogeneity 

Previous studies have reported that association with DNA 
methylation measured in whole blood samples can be driven by 
blood cell subtype heterogeneity [25]. To adjust for this, 
eosinophil, lymphocyte, neutrophil and monocyte cell counts in 
the blood samples were included (as fixed effects) in the full and 
null models for the ten most highly associated probes. 106 out of 
115 subjects had complete blood cell counts available and were 
included in this analysis. 



Exploring methylation changes in monozygotic twins 

Monozygotic twin pairs with PTA and Illumina HumanMethy- 
lation 27 k data were selected for the MZ discordance analysis 
(n = 2 1 pairs). Intra-pair DNA methylation difference per probe 
was calculated as the difference in DNA methylation residuals 
(adjusted for chip and position on the chip) between co-twins. 
DNA methylation residuals were calculated from quantile 
normalised fi values per probe. Differences in DNA methylation 
were compared to differences in PCI were using Spearman rank 
correlation. 

Effect of DNA methylation on gene expression 

To investigate the influence of DNA methylation on gene 
expression, expression levels in skin tissue collected as part of the 
Multiple Tissue Human Expression Resource (MuTHER) (http:// 
www.muther.ac.uk) were examined [26]. Quantile normalised 
gene expression in skin was adjusted for experimental batch effect 
and RNA concentration in the tissue sample and residuals 
correlated with DNA methylation residuals (adjusted for chip 
and position on the chip) at the corresponding probes using 
Pearson correlation. Furthermore, skin expression residuals were 
correlated with PC 1 values, to test for an effect of gene expression 
on the phenotype. 

Results 

Subjects and phenotypes 

Two independent samples with hearing data and DNA 
methylation profiles were selected from the TwinsUK registry to 
perform the discovery EWAS (n = 1 1 5) using Illumina Human- 
Methylation 27 k profiles, and the replication EWAS (n = 203) 
using Illumina HumanMethylation 450 k profiles. Subjects 
included in the discovery EWAS had a mean age of 56.7 years 
(±7.9 years of standard deviation from the mean, age range 33-80 
years) and included 25 dizygotic twin (DZs) pairs, 2 1 monozygotic 
twin (MZs) pairs and 23 singletons. The replication sample 
included 203 females, comprising 61 MZ twin pairs, 22 DZ twin 
pairs and 37 unpaired twins (singletons), with a mean age of 63.21 
(± 8.87 years of standard deviation from the mean, age range 41— 
82 years). The discovery and replication samples are summarised 
in Table 1. Variance in PTA was summarised using principal 
component analysis, where PCI represented the threshold shift 
over all frequencies (0.125-8.0 kHz) and captured 54.25% of the 
variance. A high PCI score thus corresponded to reduced hearing 
ability [27]. 

DNA methylation profiles 

Genome-wide DNA methylation levels were obtained in the set 
of 1 15 female twins using the Illumina 27 k array. The majority of 
autosomal CpG sites included in this analysis were unmethylated 
(P<0.3, 68.9% of probes), few probes were hemi-methylated ((5: 
0.3-0.7, 11.2% of probes) or fully methylated (P>0.7, 19.9% of 
probes). 

Discovery EWAS 

Genome-wide DNA methylation levels at 24,461 autosomal 
probes were previously obtained in the set of 1 1 5 discovery female 
twins using the Illumina 27 k array [1 1] and compared to hearing 
ability. DNA methylation at 2,519 (out of 24,641) probes was 
nominally associated (ANOVA p-value^O.05) with hearing ability 
for PC 1 . A Manhattan plot of the EWAS for hearing PC 1 is shown 
in Figure 1 and the strongest signal reached Bonferroni adjusted 
genome-wide suggestive evidence for association (p = 6.9 x 10~ 6 ). 
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The most highly associated probe was cgO 1161216 which maps 
to the promoter region of transcription factor 25 (TCF25 
)(beta±se = -0.245±0.05, p = 6.6xl0~ 6 ). Further associations 
were observed for CpG sites in the promoter regions of the 
phosphoglucomutase 3 (PGM3) gene (beta±se = — 0.26±0.06, 
p = 4.5x10 J ), the cysteine dioxygenase type 1 (CD01) gene 
(beta±se= -0.24±0.06, p = 4.7xl0~ 5 ), the nucleolar complex 
associated 2 homolog (NOC2L) gene (beta±se = -0.20±0.05, 
p = 5.4xl0~ 5 ), the myosin binding protein C (MYBPC3) gene 
(beta±se= -0.19±0.05, p = 5.4xl0~ 5 ), the fibroblast growth 
factor receptor 1 (FGFR1) gene (beta±se = -0.24±0.06, 
p = 5.7 x 10 - '), the DNA polymerase epsilon catalytic subunit 
(POLE) gene (beta±se = -0.16±0.04, p = 6.3xl0" 5 ), vacuolar 
protein sorting 4 homolog B (VPS4B) gene (beta±se = 0.20±0.05, 
p = 6.5xl0 -5 ), the heterogeneous nuclear ribonucleoprotein AO 
(HNRNPAO) gene (beta±se = 0.14±0.03, p = 6.9xl0" 5 ), and 
probe cg25017250 (beta±se = -0.23±0.06, p = 7.0xl0" 5 ) 
mapping to the apolipoprotein C-4 (AP0C4) gene. The ten most 
highly associated EWAS probes are listed in Table 2. 

After exclusion of chronological age as a fixed effect, association 
of DNA methylation with hearing PC 1 remained significant for all 
of the ten most highly associated probes (Table 2). 

Replication of highly associated EWAS probes 

The ten most highly associated CpG probes from the discovery 
sample were examined in the replication sample (Table 1). 
Association between DNA methylation and PCI was replicated 
at 2 probes - in the promoter regions of genes TCF25 and POLE 
(Table 2 and depicted in Figure 2). Figure 2 depicts the associa- 
tion between raw methylation betas with hearing PC 1 at TCF25 
and POLE in the discovery and replication samples. While probe 
cgOl 161216 (TCF25) was hypomethylated (B<0.3) in all subjects, 
probe cgl8877514 (POLE) was hypermethylated (B>0.7) 
(Figure 2). The association between adjusted DNA methylation 
residuals and PC 1 at TCF25 and POLE in the discovery and 
replication samples can be found in Figure SI. None of the 
replicating DMRs showed an underlying association of single 
nuclear polymorphisms with PC 1 200 kb up- and downstream of 
the respective genes (TCF25, POLE) in a genome-wide association 
study. 

After exclusion of chronological age as a fixed effect, association 
of DNA methylation with PC 1 remained significant at TCF25 and 
POLE (cgOl 161216: p(no age) = 1.06 x 10" 8 ; cgl8877514 p(no 
age) = 2.83 xlO" 2 )(Table 2). 

To assess the behaviour of additional probes mapping to the 
TCF25 and POLE loci with respect to hearing, the association 
between DNA methylation and hearing PC 1 was explored for all 
probes mapping to TCF25 and POLE according to hgl9 
(Figure 3). According to the locus plots, 3 further nominally 
significant associated DMRs (p<0.05) mapped to each TCF25 
and POLE. 

Although DNA methylation was not found significantly 
associated at 8 out of 10 probes in the replication sample, DNA 
methylation at five further probes (cg25383093, cgl9923810, 
cg21370143, cgl5791248 and cg25017250) showed the same 
direction of effect as in the discovery sample (Table 2). In the 
meta-analysis of results from the 27 k and 450 k chips, DNA 
methylation at 7 out of 10 probes was nominally significantly 
associated with PC 1 (Table 2), with differential DNA methylation 
at TCF25 showing the most significant association (cgO 1161216, 
p = 4.89xlO" 9 )(Table 2). 
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Figure 1 . Manhattan Plot of PCI EWAS results. The manhattan plot depicts the significance of association with PC1 as the negative logarithm 
of the p-value (-log(p-value)) versus the chromosomal location (chromosomes) for each of the 24,641 tested DNA-methylation probes. The red line 
defines a Bonferroni adjusted genome-wide suggestive significance threshold of p = 6.9x10~ 6 . The ten most highly associated probes are located 
above the horizontal blue line corresponding to p<6.985x10~ . 
doi:10.1371/journal.pone.0105729.g001 



Validation of TCF25 using MeDIPseq 

To validate the peak EWAS DMR using a different technology, 
TCF25 DNA methylation levels based on MeDIPseq data were 
also explored for association with hearing in 46 unrelated females 
from TwinsUK [24]. The mean age of subjects in the validation 
sample was 62.28 (±7.86 years of standard deviation from the 
mean, age range 43-86 years) (Table 1). MeDIPseq DNA 
methylation levels at a 1 kb locus overlapping probe 
cgO 1161216 were selected and compared to PC 1 . DNA methyl- 
ation at this locus was significantly associated with hearing PCI 
(p = 4.09x10 ) and showed the same direction of effect 
(beta±se = -8.72 x 10" B ±4.13 x 10" 6 ) as both the discovery 
EWAS and replication datasets. 

Blood cell heterogeneity 

To account for potential effects of blood cell heterogeneity, the 
peak DMRs were also explored for association with proportion of 
eosinophils, lymphocytes, neutrophils and monocytes in a subset of 
106 subjects from the discovery sample. The ten most highly 
associated probes in the discovery EWAS remained significantly 
associated (p>0.005) with PCI after adjustment for blood cell 
heterogeneity. 

Monozygotic co-twin study 

MZ discordance analyses were performed in 2 1 female MZ twin 
pairs (n = 42) selected from the discovery sample with a mean 
age of 55.43 years (±6.93 years of standard deviation, age range: 
45-68 years) (Table 1). Mean intrapair difference in PCI was 
-0.42 (±1.34 sd, range: 3.47 and -2.86). The intra-pair 



differences in PC 1 were compared with intra-pair differences in 
DNA methylation at 24,641 autosomal CpG sites. Of these, 
794 CpG sites were nominally significant (p<0.05). PCI discor- 
dance was most strongly associated with differential methylation at 
lysophosphatidic acid phosphatase 6 (ACP6, cg01377755, r= — 
0.75, p= 1.2x10 4 ). Further strongly correlated differentially 
methylated genes included myocyte enhancer factor 2D (MEF2D, 
cg08 156349, r=-0.75, p= 1.4xl0~ 4 ), tachykinin precursor 1 
{TAC1, cg07550362, r=-0.72, p = 3.3xl0~ 4 ), ATPase family 
AAA domain-containing 3C (ATAD3C, cg27383362, r=-0.70, 
p = 5.5x10 4 ), brain-specific serine protease 3 (PRSS12, 
cg21208104, r = 0.70, p = 6.3xl0" 4 ), ADAM metallopeptidase 
domain 18 (ADAM18, cg23566335, r = 0.70, p = 6.5xl0" + ), 
chromobox homolog 2 (CBX2, cg22892904, r=-0.69, 
p = 7.8xl0" 4 ), septin 3 (SEPT3, cg04283938, r=-0.68, 
p = 8.6xl0" 4 ), transmembrane protein 121 (TMEM121, 
cg23886551, r=-0.68, p = 8.6xl0" 4 ) and torsin family 1 
member B (TOR1B, cgl4299800, r=-0.68, p = 9.1xl0" 4 ) 
(Table 3). 

Influence of DNA methylation on gene expression in skin 

To investigate the influence of DNA methylation on gene 
expression, expression profiles in skin were explored because skin 
originates from the same embryonic tissues as the inner ear, and 
expression profiles were not available for the cochlea. For 172 
individuals with 27 k array data, DNA methylation at the two 
replicating probes (cgO 1161216 and eg 1887 75 14) was examined 
for association with gene expression (TCF25 and POLE, 
respectively). After adjustment of both DNA methylation and skin 
expression for batch effects, DNA methylation residuals and gene 
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PC1 versus raw DNA methylation at cg01161216 



• 27K 

• 450K 



**•« 



B 



PC1 versus raw DNA methylation at cg18877514 



• 27K 

• 450K 




Figure 2. Association of hearing PCI values and raw DNA methylation at TCF25 (cgOl 161216) and POLE (eg 1887751 4). A, B Hearing 
PC1 values were plotted versus raw DNA methylation betas for both the discovery (27 k, red dots) and the replication (450 k, blue dots) samples. 
Linear regression lines were fitted for both datasets (27 k:red line, 450 k:blue line). 
doi:1 0.1 371 /journal.pone.01 05729.g002 



expression residuals showed a weak negative correlation for 
TCF25 (r = -0.02) (Figure 4, A) and POLE (r = -0.06) (Figure 4, 
B). In general, DNA methylation in whole blood was only weakly 
correlated with gene expression in skin tissue. Furthermore, the 
effect of gene expression levels of TCF25 and POLE on hearing 
ability was explored. Gene expression showed a weak positive 
correlation with PCI values (TCF25: r = 0.12; POLE: 
r = 0.16)(Figure 4, Panel C and D). 



Discussion 

Changes in DNA methylation have been associated with 
increasing age and age-related disorders [11]. Here, for the first 
time the effect of genome-wide DNA methylation on hearing 
ability was investigated. Genome-wide association and candidate 
gene studies of hearing ability with age have yet to explain much of 
the estimated variance in this phenotype. Our approach identified 



Locus plot of DNA methylation at TCF25 



B 



Locus plot of DNA methylation at POLE 




1 1 1 1 1 — 1 1 1 1 1 1 1 

89940000 89950000 89960000 89970000 89980000 133200000 133220000 133240000 133260000 

Position on chr 16 [bp] Position on chr 12 [bp] 

Figure 3. DNA methylation association with hearing PCI at 7C;F.25and POLE. Hearing PC1 association with all available lllumina 450 k DNA 
methylation probes annotated to TCF2S (A) and POLE (B) gene regions according to hg19. Significance of association with hearing PCI expressed as 
the negative logarithm of the p-value was plotted against base pair [bp] location for chromosomes 16 and 12, respectively. The red horizontal lines 
represent the significance threshold for nominal significance (p = 0.05). 
doi:1 0.1 371 /journal.pone.01 05729.g003 
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Table 3. Results of the MZ intra-pair difference association analysis. 





MZ pair difference analysis (n = 42) 








p-value 


cgOl 377755 


ACP6 


-0.753 


1 .24E-04 


cg081 56349 


MEF2D 


-0.749 


1.41 E-04 


cg07550362 


TAC1 


-0.722 


3.25E-04 


cg27383362 


ATAD3C 


-0.703 


5.49E-04 


cg21208104 


PRSS12 


0.697 


6.26E-04 


cg23566335 


ADAM18 


0.696 


6.47 E-04 


cg22892904 


CBX2 


-0.688 


7.82E-04 


cg04283938 


SEPT3 


-0.684 


8.59E-04 


cg23886551 


TMEM121 


-0.684 


8.59E-04 


eg 14299800 


TOR1B 


-0.682 


9.1 3 E-04 



This table shows the results for the MZ discordance analysis. Results are listed for the ten most highly correlated probes with corresponding gene, Spearman rank 
correlation coefficient (rho) and significance of correlation (p-value). 
doi:1 0.1 371/joumal.pone.01 05729.t003 



epigenetic changes at a number of genes that were associated with 
hearing ability, and two of these changes in genes TCF25 and 
POLE replicated in an independent sample. DNA methylation 
levels at the strongest signal in TCF25 validated using an 
alternative method (MeDIPseq). Hearing PCI was also strongly 
associated with DNA methylation at FGFR1, a gene known to be 
essential for maintenance of glial cells and cochlear neurons in the 
spiral ganglion [28]. These findings suggest that epigenetic 
changes may account for the variance in severity and age of 
onset of ARHI. 

Nominally significant associations (p<0.05) with PCI, which 
represents the overall threshold shift in the pure-tone audiogram 
and hence impaired hearing ability, were identified at 2,519 CpG 
sites. The ten most highly associated probes remained nominally 
significant after exclusion of chronological age as a fixed effect in 
the model, showing that none of these associated probes are age- 
related differentially methylated probes. Furthermore association 
remained significant after adjustment for blood cell heterogeneity, 
indicating that blood cell subtypes were not driving these 
association signals. Two of the ten signals were replicated in a 
second independent sample. Changes in DNA methylation in the 
promoter region of TCF25 were highly associated with PCI in 
both the discovery EWAS and the replication cohort, with meta- 
analysis p = 4.89 xlO -9 . The meta-analysis of the discovery and 
replication sample was conducted to determine the joint effect of 
both samples; nevertheless the results of this analysis were driven 
primarily by the discovery EWAS findings. 

As both the discovery and replication data used the same array 
design from Illumina based on DNA hybridisation an alternative 
technique, MeDIPseq, was used to validate our findings. 
MeDIPseq in venous blood from 46 unrelated samples confirmed 
the association between hearing PC 1 and DNA methylation levels 
at TCF25 (p = 0.04). This transcription factor belongs to the 
family of basic helix-loop-helix transcription factors, which is 
widely expressed in many organs including dorsal root ganglia in 
mouse embryos [29]; however mouse models of Tcf25 deficiency 
are not yet available. Over-expression of Tcf25 leads to increased 
cell death and binding to the X-linked inhibitor of apoptosis 
protein [30]. Using the ENCODE database [31], probe 
cgO 1161216 maps to an area with enhancer and promoter 
associated histone marks and transcription factor binding sites. 
Differential expression of TCF25 might be involved in increased 



cell death of sensory cells and neurons of the cochlea, resulting in 
ARHI. 

We also identified a differentially methylated DNA methylation 
probe mapping to the promoter of the POLE gene. POLE is a 
DNA polymerase essential for elongation of the leading strand in 
cell division. In addition, POLE is involved in cell cycle regulation 
and therefore regulates a variety of cellular processes. According to 
UCSC browser and the ENCODE database [31], probe 
cgl8877514 maps to an area rich in repeating elements with 
enhancer and promoter associated histone marks and transcription 
factor binding sites. Pole knockout mice with a random gene 
disruption are embryonic lethal, while Pole targeted knock-in mice 
present with premature death due to cancer and increased 
tumourigenesis in general [32]. 

Among the top associations in the discovery EWAS was a DMR 
in the promoter of FGFR1 , a gene known to be essential for 
maintenance of glial cells and cochlear neurons in the spiral 
ganglion [28]. However, our replication study did not confirm the 
differential methylation in the promoter of FGFR1 identified in 
the discovery EWAS, but the DMR did manifest the same 
direction of effect. This gene is of particular interest, having been 
associated with hearing ability in mice [28]. FGFR1 encodes a 
fibroblast growth factor receptor, reported to be essential for 
healthy development of the organ of Corti [33]. Conditional 
knockout of fibroblast growth factor receptors (FGFR1 and 
FGFR2) in glial cells in the spiral ganglion resulted in loss of 
spiral ganglion neurons and age-related hearing loss in mice [28]. 
Our results show a negative association between DNA methylation 
at the promoter of FGFR1 and hearing PCI (beta= -0.24±0.06 
se), indicating that greater methylation (and expected reduced 
gene expression) of FGFR1 showed good hearing ability. This 
direction of effect is not consistent with that observed in mouse 
cochlea [28]. 

DNA methylation in the promoter of genes has been associated 
with repression of gene expression. At our peak DMR in TCF25, 
DNA methylation was found minimally negatively correlated with 
gene expression in skin (r= —0.02; POLE: r= —0.06) - the tissue 
with the most embryologic similarity to cochlea. However, DNA 
methylation may be highly tissue specific [34,35] and in this study 
DNA methylation was determined from whole blood samples. 
That gene expression in skin showed a weak positive correlation 
with hearing PC 1 for both TCF25 (r = 0. 1 2) and POLE (r = 0. 1 6), 
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DNA methylation at cg01 161216 versus TCF25 expression in skin 



DNA methylation at cg18877514 versus POLE expression in skin 
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Figure 4. Effect of DNA methylation on gene expression in skin and effect of gene expression on PCI. A. DNA methylation residuals 
showed a weak negative correlation (r= -0.02) with expression residuals of TCF25 in skin samples. Both quantile normalised DNA methylation betas 
and quantile normalised gene expression values were adjusted for experimental batch effects (chip and position on the chip for methylation betas 
and experimental batch and RNA concentration for gene expression profiles) previous to analysis. The regression line (blue line) depicts the linear 
association between DNA methylation residuals and gene expression residuals. B. DNA methylation residuals at probe cg18877514 were weakly 
negatively correlated (r=— 0.06) with POLE expression residuals in skin tissue. Both quantile normalised DNA methylation betas and quantile 
normalised gene expression values were adjusted for experimental batch effects (chip and position on the chip for methylation betas and 
experimental batch and RNA concentration for gene expression profiles) prior to analysis. The regression line (blue line) depicts the linear association 
between DNA methylation residuals and gene expression residuals. C. TCF25 expression residuals in skin showed a weak positive correlation (r = 0.1 2) 
with PC1 . Quantile normalised gene expression values were adjusted for experimental batch effects and RNA concentration. The regression line (blue 
line) depicts the linear association between gene expression residuals and PCI values. D. POLE expression residuals in skin showed a weak positive 
correlation (r = 0.1 6) with PCI . Quantile normalised gene expression values were adjusted for experimental batch effects and RNA concentration. The 
regression line (blue line) depicts the linear association between gene expression residuals and PCI values. 
doi:1 0.1 371 /journal.pone.01 05729.g004 



indicates that individuals with decreased hearing ability (high PC 1 
value) show higher RNA levels of TCF25 and POLE in skin. 
Whether these findings pertain to RNA expression in the inner ear 
remains to be determined. 



Monozygotic twin pairs are a preferred study sample for 
epigenetic studies as they are assumed to be genetically identical. 
In addition, both dizygotic and monozygotic twin pairs show an 
increased proportion of shared environment due to the nature of 
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their time shared in uterus and upbringing. The MZ discordance 
analysis was performed to best utilise the unique study sample 
presented here and for completeness. Nevertheless, the relatively 
low sample size and restricted discordance within the twin pairs 
limited the statistical power to detect strong epigenetic effects. 
Association was examined between intra-pair discordance for PC 1 
and intra-pair differences in DNA methylation at CpG sites 
genome-wide. The most highly associated probes were found in 
the promoters of ACP6 and MEF2D. The function of acid 
phosphatase 6 is yet unknown and Acp6 knockout mice are 
described as phenotypically normal [32]. In contrast, myocyte 
enhancer factor 2D is a member of the myocyte enhancer factor 
family of transcription factors, which are involved in neuronal 
development and differentiation under regulation of class 2 histone 
deacetylases. MEF2D is expressed in mouse cochlear neurons and 
sensory cells at P15 and was diminished in IGF knockout mice, 
which show sensorineural hearing loss [36]. The data indicate that 
MEF2D is a plausible candidate gene for ARHI and may be 
under epigenetic control. 

Our study has several strengths and limitations. DNA methyl- 
ation is likely to play an important role in gene expression 
contributing to important phenotypic differences between tissues, 
between individuals and with age. Methods of analysis of 
methylation data are in their infancy: there are many important 
covariates to be considered. We elected to remove one of these, 
gender, by confining our studies to females, which predominate in 
the TwinsUK database. Thus our results pertain to women and 
may not extrapolate to men. Strengths included ability to exclude 
age and blood cell heterogeneity as potential confounders. The 
high proportion of related individuals in this sample reduced both 
the genetic and environmental variance compared to a population 
sample of unrelated individuals. Although the discovery and 
replication datasets were well matched for gender, ethnicity, age 
and hearing ability, the replication sample included by chance a 
higher proportion of monozygotic twin pairs (450 k sample: 60% 
MZs) compared to the discovery sample (27 k sample: 37% MZs), 
which might have resulted in the reduced significance of 
associations obtained in the replication sample. Association in 
the EWAS did not reach epigenome-wide significance by 
Bonferroni corrected significance levels (considering 24,641 
independent tests: p£2.03xl0 6 ). However, DNA methylation 
of neighbouring CpG sites is unlikely to be independent thus a 
Bonferroni correction may be considered overly stringent. Taking 
co-methylation into account by correcting for the number of 
genes, the peak DMR in TCF25 surpassed genome-wide 
suggestive evidence for association. Further limitations of the 
study included the choice of tissue: although DNA methylation is 
tissue specific, whole blood samples were used as an initial 
approach to this investigation because they were readily available 
and inner ear tissue from humans was not. In addition, 
discordance in hearing ability within TwinsUK monozygotic twin 
pairs was relatively limited. Finally, it should be noted that this 
study makes no assumptions about causal relationships between 
DNA methylation and ARHI. A longitudinal study design would 
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be required to confirm that the methylation changes inferred by 
these results predated the onset of hearing impairment. 

In conclusion, this is the first study investigating the association 
between hearing ability with age and DNA methylation genome- 
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promoters of 1 0 genes were identified, of which two (TCF25 and 
POLE) were replicated in an independent cohort. Functional 
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in the development of ARHI would highlight the impact of 
changes in DNA methylation with age and therefore be of 
fundamental importance not only for hearing loss but also other 
age-related disorders. 
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replication (450 k, blue dots) samples. A linear regression lines was 
fitted for both datasets (27 k:red line, 450 k:blue line). 
(TIF) 

Acknowledgments 

This work was funded by Action on Hearing Loss and AgeUK. 

The authors would like to thank all volunteers from the TwinsUK 
register, who contributed to this study and acknowledge the work of all 
researchers who helped in the completion of this work. The study was 
funded by the Wellcome Trust; European Community's Seventh 
Framework Programme (FP7/2007-2013). The study also receives support 
from the National Institute for Health Research (NIHR) BioResource 
Clinical Research Facility and Biomedical Research Centre based at Guy's 
and St Thomas' NHS Foundation Trust and King's College London. Tim 
Spector is holder of an ERC Advanced Principal Investigator award. 

Author Contributions 

Conceived and designed the experiments: LEW JTB FMKW PT CJS. 
Performed the experiments: LEW PT CJS. Analyzed the data: LEW. 
Contributed reagents/materials/analysis tools: PD TDS CJS. Wrote the 
paper: LEW. 



5. Rakyan VK, Down TA, Balding DJ, Beck S (201 1) Epigenome-widc association 
studies for common human diseases. Nat Rev Genet 12: 529-541. 

6. Okano M, Bell DW, Haber DA, Li E (1999) DNA Methyltransferases Dnmt3a 
and Dnmt3b Are Essential for De Novo Methylation and Mammalian 
Development Cell 99: 247-257. 

7. Bestor I'H (2000) The DNA methyltransferases of mammals. Human Molecular 
Genetics 9: 2395-2402. 

8. Wong CCY, Caspi A, Williams B, Craig IW, Houts R, ct al. (2010) A 
longitudinal study of epigenetic variation in twins. Epigcnetics 5: 516-526. 



PLOS ONE | www.plosone.org 



10 



September 2014 | Volume 9 | Issue 9 | e1 05729 



Epigenome-Wide DNA Methylation in Hearing Ability 



9. Bell JT, Pai AA, Pickrcll JK, Gaffney DJ, Piquc-Rcgi R, ct al. (2011) DNA 
methylation patterns associate with genetic and gene expression variation in 
HapMap cell lines. Genome Biol 12: RIO. 

10. BellJT, Spector TD (201 1) A twin approach to unraveling epigenetics. Trends in 
Genetics 27: 116-125. 

11. BellJT, Tsai P-C, Yang T-P, Pidslcy R, NisbctJ, et al. (2012) Epigenome-Wide 
Scans Identify Differentially Methylated Regions for Age and Age-Related 
Phenotypes in a Healthy Ageing Population. PLoS Genet 8: el002629. 

12. Dempster EL, Pidsley R, Schalkwyk LC, Owens S, Georgiades A, et al. (201 1) 
Disease-associated epigenetic changes in monozygotic twins discordant for 
schizophrenia and bipolar disorder. Human Molecular Genetics 20: 4786-4796. 

13. Rakyan VK, Down TA, Maslau S, Andrew T, Yang T-P, et al. (2010) Human 
aging-associated DNA hypcrmcthylation occurs preferentially at bivalent 
chromatin domains. Genome Research 20: 434-439. 

14. Gruickshanks KJ, Wiley XL, Xweed XS, Klein BE, Klein R, et al. (1998) 
Prevalence of hearing loss in older adults in Beaver Dam, Wisconsin. Xhe 
Epidemiology of Hearing Loss Study. Am J Epidemiol 148: 879-886. 

15. Provenzano MJ, Domann FE (2007) A role for epigenetics in hearing: 
Establishment and maintenance of auditory specific gene expression patterns. 
Hearing Research 233: 1-13. 

16. Friedman L, Avraham K (2009) MicroRNAs and epigenetic regulation in the 
mammalian inner ear: implications for deafness. Mammalian Genome 20: 581— 
603. 

17. Wilkin DJ, Libcrfarb R, Davis J, Levy HP, Cole WG, et al. (2000) Rapid 
determination of COL2A1 mutations in individuals with Stickler syndrome: 
analysis of potential premature termination codons. American Journal of 
Medical Genetics 94: 141-148. 

18. Donoso LA, Edwards AO, Frost AX, Ritter R, Ahmad N, ct al. (2003) Clinical 
variability of Stickler syndrome: role of exon 2 of the collagen GOL2A1 gene. 
Survey of ophthalmology 48: 191-203. 

19. Buschdorf JP, Strading WH (2004) A WW domain binding region in methyl- 
CpG-binding protein MeCP2: impact on Rett syndrome. Journal of Molecular 
Medicine 82: 135-143. 

20. BSA (201 1) Rccommende Procedure Pure-tone air-conduction and bone- 
conduction threshold audiometry with and without masking. Reading: British 
Society of audiology. 

21. Bibikova M, LeJ, Barnes B, Saedinia-Melnyk S, Zhou L, et al. (2009) Genome- 
wide DNA methylation profiling using lnfinium assay. Epigenomics 1 : 1 77-200. 

22. Bibikova M, Barnes B, Xsan C, Ho V, Klotzle B, et al. (201 1) High density DNA 
methylation array with single CpG site resolution. Genomics 98: 288-295. 

23. Wilier CJ, Li Y, Abecasis GR (2010) MEXAL: fast and efficient meta-analysis of 
genomewide association scans. Bioinformatics 26: 2190-2191. 



24. Bell JX, Loomis AK, Butcher LM, Gao F, Zhang B, et al. (2014) Differential 
methylation of the XRPA1 promoter in pain sensitivity. Nat Commun 5. 

25. Horvath S, Zhang Y, Langfeldcr P, Kahn RS, Boks MP, et al. (2012) Aging 
effects on DNA methylation modules in human brain and blood tissue. Genome 
Biol 13: R97. 

26. Grundberg E, Small KS, Hedman AK, Nica AG, Buil A, et al. (2012) Mapping 
cis- and trans -re gulatory effects across multiple tissues in twins. Nat Genet 44: 
1084-1089. 

27. Huyghe JR, Van Laer L, HendrickxJJ, Fransen E, Demecstcr K, et al. (2008) 
Genome-wide SNP-based linkage scan identifies a locus on 8q24 for an age- 
related hearing impairment trait. Am J Hum Genet 83: 401—407. 

28. Wang SJ, Furusho M, D'Sa G, Kuwada S, Conti L, et al. (2009) Inactivation of 
fibroblast growth factor receptor signaling in myelinating glial cells results in 
significant loss of adult spiral ganglion neurons accompanied by age-related 
hearing impairment. Journal of Neuroscicnce Research 87: 3428—3437. 

29. Olsson M, Durbeej M, Ekblom P, Hjalt X (2002) Nulpl, a novel basic helix- 
loop-hclix protein expressed broadly during early embryonic organogenesis and 
prominently in developing dorsal root ganglia. Cell and Xissue Research 308: 
361-370. 

30. Steen H, Lindholm D (2008) Nuclear localized protein-1 (Nulpl) increases cell 
death of human osteosarcoma cells and binds the X-linked inhibitor of apoptosis 
protein. Biochemical and Biophysical Research Communications 366: 432-437. 

31. Rosenbloom KR, Sloan CA, Malladi VS, Dreszer XR, Learned K, et al. (2013) 
ENCODE Data in the UCSC Genome Browser: year 5 update. Nucleic Acids 
Research 41: D56-D63. 

32. EppigJX, Blake JA, Bult CJ, Kadin JA, Richardson JE (2012) Xhe Mouse 
Genome Database (MGD): comprehensive resource for genetics and genomics 
of the laboratory mouse. Nucleic Acids Research 40: D881-D886. 

33. Pirvola U, Ylikoski J, Xrokovic R, Hebert JM, McConnell SK, et al. (2002) 
FGFR1 Is Required for the Development of the Auditory Sensory Epithelium. 
Neuron 35: 671-680. 

34. Byun H-M, Siegmund KD, Pan F, Weisenberger DJ, Kanel G, et al. (2009) 
Epigenetic profiling of somatic tissues from human autopsy specimens identifies 
tissue-and individual-specific DNA methylation patterns. Human Molecular 
Genetics 18: 4808-4817. 

35. Ladd-Acosta C, Pevsner J, Sabunciyan S, Yolken RH, Webster MJ, et al. (2007) 
DNA methylation signatures within the human brain. Xhe American Journal of 
Human Genetics 81: 1304-1315. 

36. Sanchcz-Caldcron H, Rodriguez-dc La Rosa L, Milo M, Pichel JG, Holley M, 
et al. (2010) RNA microarray analysis in prenatal mouse cochlea reveals novel 
IGF-I target genes: implication of MEF2 and FOXM1 transcription factors. 
PLoS One 5: c8699. 



PLOS ONE | www.plosone.org 



11 



September 2014 | Volume 9 | Issue 9 | e105729 



